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We study the dependence of the transport properties of square Josephson Junctions arrays 
with the direction of the applied dc current, both experimentally and numerically. We present 
computational simulations of current- voltage curves at finite temperatures for a single vortex in the 
array (Ha /$o = / = 1/L ), and experimental measurements in 100 x 1000 arrays under a low 
magnetic field corresponding to / ps 0.02. We find that the transverse voltage vanishes only in the 
directions of maximum symmetry of the square lattice: the [10] and [01] direction (parallel bias) and 
the [11] direction (diagonal bias). For orientations different than the symmetry directions, we find a 
finite transverse voltage which depends strongly on the angle <j> of the current. We find that vortex 
motion is pinned in the [10] direction (<f> = 0), meaning that the voltage response is insensitive to 
small changes in the orientation of the current near <f> — 0. We call this phenomenon orientational 
pinning. This leads to a finite transverse critical current for a bias at cj> — and to a transverse 
voltage for a bias at <j> ^ 0. On the other hand, for diagonal bias in the [11] direction the behavior is 
highly unstable against small variations of 4>, leading to a rapid change from zero transverse voltage 
to a large transverse voltage within a few degrees. This last behavior is in good agreement with our 
measurements in arrays with a quasi-diagonal current drive. 



I. INTRODUCTION 



The interaction between the periodicity of vortex lattices, j^VX) and periodic pinning potentials in pinerconduc- 
tors has raised a great interest both in equilibrium systemsciij and in driven non-equilibrium systems .EjLj Several 
techniques have been used to artificially fabricate periodic pinning in superconducting samples: thickness modulated 
filmsjj wire networks,0 Josepison junction arrays,0'u magnetic dot arraysjj sub-micron hole latticeso and pinning 
induced by Bitter decoration.D 

The ground states of these systems, which result from the competition between the vortex-vortex and the vortex- 
pinning|-|iateractions, can be either commensurate or incommensurate vortex structures depending on the vortex 
density.EJiJ These conmensurability effects in the ground state vortex configurations lead to enhanced critical currents 
and resistance minima for the "matching" and for the "fractional" (submatching) vortex densities where the VL is 
strongly pinned. At finite temperatures, it has been shown that there are both a depinuing and a melting transition, 
which can occur either sequentially or simultaneously depending on the magnetic field.lj 

Very recently- the physics of driven vortices .under periodic pinning has been studied numerically both at zero 
temperaturetj'Eil and at finite temperatures .E3tj At T — there is a complex variety of dynamic phases lij At finite 
T there are two dynamic transitions when increasing temperature at high drives: there is first a transverse depinning 
and second a melting transition of the moving vortex latticeEj 

Most of the effects of periodic pinning that have been studied are related to conmensurability phenomena and 
the breaking of translational symmetry in these systems. Less studied is the effect of the breaking of rotational 
symmetry in periodic pinning potentials, in particular regarding transport properties. One question of interest is how 
the motion of vortices changes when the direction of the driving current is varied, ff there is rotational symmetry, the 
vortex motion and voltage response should be insensitive to the choice of the direction of the current. However, it is 
clear that in a periodic pinning potential the dynamics may depend on the direction of the current. For example, in 
square Josephson junction arrays (JJA) it has been found that the existence of fractional giant Shapiro steps (FGGS) 
depends on the orientation of the current bias. When the JJA isjdriven in the [11] direction the FGGS are absent, 
while they are very large when the drive is in the [10] directionO Another example of more recent interest is the 
phenomenon of transverse critical current in superconductors with pinning.lij It has been found that for a VL a driven 
with a high, ettrrent there is a transverse critical current when an additional small bias is applied in the perpendicular 
direction.Ejtj Furthermore, whence transverse bias is increased it is possible to have a rich behavior with a Devil's 
staircase in the transverse voltageJlil 



1 



In this paper we will study in detail the breaking of rotational invariance in square JJA. In this case, the discrete 
lattice of Josephson junctions induces a periodic egg-carton potential for the motion of vortices.u We will study here 
how the voltage response depends on the angle of the current with respect to the lattice directions of the square 
JJA. We will show that there are preferred directions for vortex motion for which there is orientational pinning. 
This leads to an anomalous transverse voltage when vortices are driven in directions different from the symmetry 
directions. An analogous effect of a transverse voltage due to the guided motion of vortices has been observed in 
YBCO superconductors with twin boundaries.LS Another related case is the intrinsic breaking of rotational symmetry 
of d-wave superconductivity which causes an angle-dependent transverse voltage for large currentsJHI Here we will 
show the differences and similarities of the square JJA with these problems. 

The paper is organized as follows. In Sec. II we present the model equations for the dynamics of the JJA, which 
will be solved in the numerical simulations. In Sec. Ill we describe the experimental details of the JJA used in the 
measurements. In Sec. IV we will present our results for both the simulations and for the experiments. In particular, 
we will present experiments corresponding to the orientation for which the effect of a transverse voltage is maximum. 
Finally in Sec. V we will compare our results with other similar effects and discuss future directions of study. 

II. MODEL 

We study—the dynamics of JJA using the resistively shunted junction (RSJ) model for the junctions of the square 
network. E3~E!I In this case, the current flowing in the junction between two superconducting islands in a JJA is modeled 
as the sum of the Josephson supercurrent and the normal current: 

^(n)=/osin^(n) + 27r ^y (1) 

where Iq is the critical current of the junction between the sites n and n + /i in a square lattice [n = (n x , n y ), fi = x, y], 
Rn is the normal state resistance and 

M (n) = 0(n + M ) - 0(n) - A^(n) = A M 0(n) - A^n) (2) 

is the gauge invariant phase difference with 
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The thermal noise fluctuations rj^ have correlations 



IkT 

(^(n,t) V (n',0) = -^<W<W5(t-t') (4) 

In the presence of an external magnetic field H we have 

A M x A M (n) = A x (n) - A x (n + y) + A y (n + x) - A y (n) 

= 2tt/, (5) 

/ = Ha 2 /§o and a is the array lattice spacing. We take periodic boundary conditions (p.b.c) in both directions in 
the presence of an external current I = (I x , I y ) in arrays with L x L junctions.Ea The vector potential is taken as 

A lt (n,t) = A° tl (n)-a ll (t) (6) 

where in the Landau gauge ^(n) = — 2tt fn y , A y (n) = and a^t) allows for total voltage fluctuations under periodic 
boundary conditions. In this gauge the p.b.c. for the phases are£j£3 

9(n x + L, n y ) = 9(n x ,n y ) 

0(n x ,n y + L) = 9(n x ,n y ) - 2tt fLn x . (7) 
The condition of a total current flowing in the x and y directions: 
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Ix - L 2 



Iy ~ L 2 



^2 !q sin 9x (n)+r) x (n, t) 



Io sin By (n) + T] y (n, t) 



h da x 
2eR N ~dT ' 



(8) 



h da 



v 



2eR N dt ' 



determines the dynamics of a p (i).@ We also consider local conservation of current, 

A„.J p (n) = ^;/ /t (n)-/ / ,(n-|i) = 0. 

After Eqs. (1,8,9) we obtain the following set of dynamical equations for the phases 



A 



2 d8(n) 



dt 
dt 



-A M -[^(n)+J7 M (n,t)] 
^-^2l][^( n ) + ^( n : *)] 



where 



S tl (n)= S m[A fl B(n)-A°Jn)-a fl ], 



(9) 

(10) 
(11) 

(12) 



we have normalized currents by Iq, time by tj = 2itcRnIo/&o, temperature by Io^o/2irkB, and the discrete Laplacian 
is 



A£0(n) = 0(n + x) + B(n - x) + 6>(n + y) + 0(n - y) 
-40(n). 



(13) 



The Langevin dynamical equations (10-11) are solved with a second order Runge-Kutta-Helfand-Greenside algo- 
rithm with time step At = 0.1 77 and integration time 5000tj after a transient of 2000rj. The discrete Laplacian is 
inverted with a fast Fourier + tridiagonalization algorithm as in Rcf. ^l|. We calculate the time average of the total 
voltage as 



V x = (v x {t)} = (da x {t)/dt) 
V y = (v y (t)) = (da y (t)/dt) 



(14) 



with voltages normalized by RnIq- 

We study the JJA under a magnetic field corresponding to a single vortex in the array, / = l/L 2 , and system sizes 
of L x L junctions, with L = 32,64. We apply a current I at an angle </> with respect to the [10] lattice direction, 



I x = I cos (j) 
Iy = I sin <j). 

Wc define the longitudinal voltage as the voltage in the direction of the applied current, 

Vi — V x cos <p + V y sin 0, 

and the transverse voltage 

Vt = —V x sin 4> + V y cos <j). 
From the voltage response, we define the transverse angle as 

tan.0t = V t /Vi 

and the voltage angle as 

tan0„ = Vy/Vx, 



(15) 
(16) 
(17) 
(18) 
(19) 



i.e. B t = B v — 4>, see Fig. 1. 

When the vortices move in the direction perpendicular to the current, there is no transverse voltage, therefore 
t = and B v = (j). 
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III. EXPERIMENTAL SETUP 



We measured Current -Voltage (IV) characteristics of square proximity-effect Pb/Cu/Pb Josephson Arrays with 
the current applied in different directions. 

The samples consist of 2500 A-thick cross-shaped Pb islands on top of a continuous 2500 A-thick copper film. 
Copper, and subsequently lead, were thermally evaporated onto a silicon substrate within the same evaporator. An 
array of 1000 x 1000 lead islands were defined by photolithographic patterning followed by an Ar ion etching. The 
cell parameter of the resulting array was 10 /xm with junctions 2 /im wide and a separation of 1 fim. A second 
photolithography step was used to define a 1 x 10 mm 2 strip with current and voltage (longitudinal and transversal) 
contacts. This mask was manually aligned in the [10] and [11] directions for different samples. 

The six-terminal measurements were made using a programable dc current source and a two channel nanovoltmctcr 
(HP 34420A) . Each channel measuring the longitudinal and transversal voltage respectively. The arrays were coooled 
up to 1.25 K in a pumped 4 He cryostat shielded by /x-metal. A superconducting solenoid was used to null the 
remaining ambient magnetic field (~ 15 mG), and apply fields to the sample. The typical periodic-in-field response 
in the resistance were observed extending in a large number of periods, which was used to determine the frustration 
applied to the sample. 

IV. NUMERICAL AND EXPERIMENTAL RESULTS 
A. Breaking of rotational invariance. 

The square lattice has two directions of maximum symmetry: the [10] and the [11] directions (and the ones obtained 
from them by 7r/2 rotations), which correspond to the directions of reflection symmetry. When the current bias is in 
the [10] direction, the angle of the current is <fi = 0, and we call it a "parallel" bias. When the current bias is in the 
[11] direction, the angle of the current is (f> = ir/4 = 45°, and we call it a "diagonal" bias. 

In the case of the parallel bias we find that the transverse voltage is zero (in agreement with the reflection symmetry) . 
This corresponds to vortex motion in the direction perpendicular to the current (0t =0). In the IV curve for the 
longitudinal voltage we find a critical current corresponding to JJpja-single vortex depinning /I 10 ' = 0.1,0 as shown in 
Fig. 2(a). Above I c it is possible to distinguish three regimes :l!irE3 (A) a single vortex regime for 0.1 < I < 0.85, 
where the IV is quasi-linear and dissipation is caused by vortex motion; (B) an intermediate crossover regime for 
0.85 < I < 1.0 and (C) a resistive regime for I > 1.0 where dissipation is caused by the Ohmic shunt resistance and 
the IV curve is linear. The vortex regime A can be described by the dynamics of a single collective degree of freedom: 
an overdamped particle moving in the periodic "egg-carton" potential,o£3 although with a non-linear viscosity.EHl The 
resistive regime C is also very simple, it can be represented by the behavior of a single junction at large currents, 
9fj,(n,t) w 2el j^ N Ifjt + 5 fJ ,(n,t). The crossover regime B is characterized by a complex collective dynamics with an 
interplay of the vortex degree of freedom with spin waves excitations, and at finite temperatures there is also a steep 
increase of vortex-antivortex excitations in this regime. 

In the case of the diagonal bias, we obtain similar results as in the parallel bias case. The transverse voltage is zero 
and therefore the vortex moves perpendicular to the current (fit = 0). The IV curve for Vi has a critical current of 
li U] = V2li W] = 0.1414. The onset of the resistive regime C is also multiplied by a factor of a/2, while the crossover 
regime B starts at nearly the same current as in the parallel bias case, see Fig. 2(b). 

For orientations different than the symmetry directions, we always find a finite transverse voltage. In order to see 
this, we study the voltage response when varying the orientation <fi of the drive while keeping fixed the amplitude / of 
the current. In Fig. 3(a) we plot the transverse angle 8 t — arctan(Vt/V/) as a function of the angle of the current cj>. 
We find that B t vanishes only in the maximum symmetry directions corresponding to angles = 0, ±45°, ±90°, . . ., 
as discussed before. Furthermore, we see that for orientations near = 0, the transverse angle basically follows 
the current angle: 9t ~ —(p. This is an indication that vortex motion is pinned in the lattice direction [10], since 
V y « 0, meaning that the voltage angle is 6 V « 0. Whenever the voltage response is insensitive to small changes in the 
orientation of the current, we will call this phenomenon orientational pinning. On the other hand, near (f> = 45° the 
transverse angle changes rapidly. We show in Fig. 3(b-c) the behavior of Q t for different currents and temperatures. 
The 9 t vs. <fi curves become smoother around <fi = 45° for increasing current as well as for increasing T. At the same 
time, the magnitude of the transverse angle decreases when increasing the current for / 3> 1, or when increasing T. 

A more direct evidence of the breaking of rotational symmetry can be seen in the parametric curves of V y (cj)) vs. 
V x (</>). In Fig. 4 we plot the values of the voltages V y and V x when varying the orientational angle 4> for different 
values of the current amplitude / and the temperature. In the case of rotational symmetry we should have a perfect 
circle. In the set of plots of Fig.4(a-c), the current amplitude is fixed and the temperature is varied. In Fig. 4(a) we 
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have I = 0.2, near the onset of single vortex motion in the regime A. In this case most of the points are either on 
the axis V x = or on the axis V y — 0, indicating strong orientational pinning in the lattice directions [10] or [01]. 
When increasing T the orientational pinning decreases and the length of the "horns" in the x and y axis decreases. 
Fig. 4(b) corresponds to I = 0.6 , near the end of regime A when the vortex is moving fast. In this case the horns 
have disappeared and orientational pinning is lost. However, the breaking of rotational symmetry is still present in 
the star-shaped curves that we find at low T. The dip at 45° in the stars are because in this direction the voltages 
are minimum, since the critical current is maximum in this case, = O.l-y/2. When increasing the temperature, 
the stars tend to the circular shape of rotational invariance. The parametric curves in the crossover regime B also 
have star shaped behavior at low T which tends to circles when increasing T. Above the onset of the resistive regime 
C the "horned" curves reappear [Fig. 4(c)]. In this case the orientational pinning corresponds to the locking of ohmic 
dissipation in the junctions in one of the lattice directions, either [10] or [01]. Once again, when increasing T the 
horned structure shrinks, and the curves evolve continuously from square shapes to circular shapes. 

The variation with current of the rotational parametric curves for a fixed temperature is shown in Figs.4(d-f). At 
a low temperature, T = 0.05, we clearly see the horned structure of the curves for almost all the currents and even 
for large currents the circular curves have "horns", see Fig. 4(d). At an intermediate temperature, T = 0.2, there are 
still some signatures of the orientational pinning [Fig. 4(e)], while for T = 1.2 all the curves are smooth and rounded 
with a slightly square shape [Fig. 4(f)]. 

B. Orientational pinning near the [10] direction. 

The orientational pinning characterizes the breaking of rotational symmetry in square arrays at low temperatures, as 
we have seen in the previous Section. When the JJA is driven in any of the "parallel" directions, <fi = 0°, 90°, 180°, 270° 
both vortex motion and dissipation are pinned along these directions. When the current is rotated a small angle, 
for example from the [10] direction, the dissipation remains pinned along the [10] direction (14^0, V y = 0). This 
effect causes a finite transverse voltage, measured with respect to the direction of the current, V t s» — V x sin(f> and a 
transverse angle 9 t = —<f>. The orientational pinning is lost at a given critical angle <p c which depends on temperature 
and current. As we can see in Fig. 3(a), for very low T and for certain values of the current, the critical angle can 
reach values very close to 45°. This leads to a large transverse voltage in arrays driven near the [11] direction as we 
will discuss later in Sec. IV. C. 

Let us analyze now the behavior of the critical angle for orientational depinning, (f> c . This can be studied by looking 
at the angle of the voltage with respect to the [10] direction, 9 V = arctan(V y /14). For (f> < (f> c we have orientational 
pinning and therefore 9 V — 0, while for <fi > 4>c we have 6 V ^0. Therefore, the onset of a finite 6 V defines the critical 
angle cj> c . In Fig. 5 we plot 9 V as a function of cj> for different currents and temperatures. There are two cases of 
interest: for a current in the single vortex regime A [Fig. 5(a)] and for a current in the resistive regime C [Fig. 5(b)]. 
We find that in both cases there is clearly a finite <j) c which decreases with temperature. In the case of the single 
vortex regime A—we see that <f> c tends to vanish at a crossover temperature corresponding to the energy scale for 
vortex depinning,Q T p i n ~ AE p i n w 0.2. The notion of "transverse critical pirrent" l c ± r has been introduced recently 
by Giamarchi and Le Doussal for driven vortex lattices in random pinning S3 When a vortex structure is moving fast, 
it can still be pinned in the transverse direction. After applying a current in the direction perpendicular to the drive, 
a finite transverse critical current-may exist at T — 0. In the case of periodic pinning, a finite I c ^ r was also proposed 
due to conmensurability effects.E-J In particular, in the periodic "egg-carton" pinning of Josephson junction arrays 
it was found in Ref. O that I c tr is finite in a wide range of temperatures for a system driven in the [10] direction. 
Here we see that due to the orientational pinning, a transverse critical current is finite only when the JJA is driven 
either in the [10] or [01] directions. In any other case it will be zero. Moreover, the critical angle in the regime A can 
be interpreted as corresponding to a transverse critical current I c ^ r — Ism(f> c for a vortex driven by a longitudinal 
current Ii = I cos <pc ■ 

In Fig. 5(b) we see that for a current in the resistive regime C there is also a critical angle for the onset of transverse 
ohmic dissipation. In this case <f> c tends to vanish at a higher temperature above the Kostcrlitz-Thouless transition, 
T KT w 0.9. 

The other lattice symmetry direction of interest corresponds to the case of "diagonal" bias, i.e. the [11] directions 
of </> = 45°, 135°, 225°, 315°. In this case, we have shown previously that the transverse voltage also vanishes. In 
agreement with this, we see both in Fig. 5(a) and Fig. 5(b) that all the curves cross in the point (<fi,9 v ) = (45°, 45°), 
which corresponds to 9 t = 6 V — cj) = 0. However, there is no orientational pinning in the [11] direction. On the 
contrary, in this direction any small deviation in the orientation of the current can cause a fast increase of the 
transverse voltage. In order to show this effect, we study the "diagonal voltage angle", which is measured with 
respect to 45°: 9' v — arctan(V^ — V x )/(V y + V x ) = 9 V — 7r/4, as a function of the deviation from the [11] direction, 



5 



Acj) — 4> — 7r/4, see Fig. 6. We find that, in contrast with the case of Fig. 5, a small deviation from the symmetry 
direction leads to a large change in the voltage angle for any current and temperature. This shows that the [11] 
direction is highly unstable against small changes in the orientation of the current, leading to an anomalously large 
transverse voltage. 

C. Transverse voltage near the [11] direction. 

In Fig. 7(a) we show our experimental voltage-current characteristics for an array of 100 x 1000 junctions at a low 
temperature T = 1.25K and at a low magnetic field. The current is applied nominally in the [11] direction, but a 
small misalignment is possible in the setup of electrical contacts, therefore <j> = 45° ± 5°. We see that for low currents 
there is a very large value of the transverse voltage Vt, which is nearly of the same magnitude as the longitudinal 
voltage V\. The transverse voltage is maximum at a characteristic current I m . Above I m , Vt decreases with increasing 
current while V\ increases. It is remarkable that these results are very different from the IV curve of Fig. 2, where 
Vt = at <f> = 45°. However, if we assume a misalignment of a few degrees with respect to the [11] direction we 
can reproduce the experimental results. In Fig. 7(b) we show the IV curves obtained numerically for <j> = 40° and 
T = 0.02. We see that for low currents Vt is close to Vt. Vt < VJ, similar to the experiment, and later Vt has a 
maximum at a current I m w 1/ cos <j> « 1.3. This corresponds to the current for which the junctions in the x-direction 
become critical {I x = 1). The range of currents we can measure experimentally is limited to the regimes B and C, 
since we can not fully access the regime A of single vortex motion due to the small voltages involved in this case. Of 
course, in the simulations we can study the full range of currents, which is shown in Fig. 7(c). Here we see that near 
the vortex depinning current the transverse voltage is also very close to Vi in a small range of currents, then when 
increasing / they separate first inside the regime A, and later in the regime B the transverse voltage approaches the 
longitudinal voltage again. 

As we saw in Fig. 3 the highest transverse voltage can be obtained for orientations near <f> — 45°. Therefore a 
slight misalignment of the array from the [11] direction is useful for studying both experimentally and numerically 
the behavior of the transverse voltage as a function of current and temperature. 

In Fig. 8 we show the dependence of 9t with current for different temperatures. The experimental results are shown 
in Fig. 8(a) where we find that 6 t first increases with current, it reaches a maximum value 9 r t nax (T) < 45° and then 
for large currents 9 t tends to zero. In Fig. 8(b) we show that the simulations with the RSJ model with <f> = 40° 
reproduce this behavior. Here we see that the maximum of the transverse angle is reached inside the regime B well 
before the onset of the resistive regime. We also find that, when going deep into the regime C, 8 t decreases with 
current: 9 t — > for I I m . In Fig. 8(c) we show that a similar behavior is obtained for other values of 4> close to 45°. 
We also observe here the full range of currents. We see that at the critical current the transverse angle 8t first has a 
maximum, then it decreases rapidly in a small range of currents, after which, for most of the regime A, 9 t is increasing 
with I before reaching a second maximum value in the crossover regime B. This shows that there are two regimes 
where the effect of anomalous transverse voltage is maximum: near the vortex depinning current, due to orientational 
pinning of vortex motion; and near the Josephson junction critical current due to the orientational "pinning" of ohmic 
dissipation. Regretfully, we can not measure the small voltages of the low current regime, therefore we were not able 
to observe experimentally the first maximum of 9 t . 

In Fig. 9 we analyze the behavior of the transverse angle as a function of temperature. We plot the value of 9 t for a 
current near the maximum value of the transverse voltage at low T. We observe experimentally that 9 t decreases with 
temperature and in particular it has a sharp decrease at T w 1.5K, as we show in Fig. 9(a). On the other hand, in the 
simulation results we see a smooth decrease of 9t with temperature [Fig. 9(b)]. The transverse angle becomes small 
at the depinning crossover temperature T p i n sa 0.2. The fact that our experimental results show a sharper decrease 
with temperature than the simulations is possibly due to vortex collective effects. The simulation results presented 
here are focused in the motion of a single vortex in the periodic pinning oLa square JJA. The vortex collective effects, 
which have to be studied for fields / > l/L 2 , will be discussed elsewhererJ 

V. DISCUSSION 

In the egg-carton potential of a square JJA there are pinning barriers for vortex motion in all the directions. The 
direction with the lowest pinning barrier is the [10] direction. Therefore the strong orientational pinning we find here 
is in the direction of the lowest pinning for motion, i.e. the direction of easy flow for vortices. The presence of a 
strong orientational pinning leads to a large transverse voltage when the systems is driven away from the favorable 
direction, to the existence of a critical angle and to a transverse critical current. On the other hand, the [11] direction 
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is the direction of the largest barrier for vortex motion in the egg-carton potential. In this case, the behavior is highly 
unstable against small variations in the angle of the drive, leading to a rapid change from zero transverse voltage to 
a large transverse voltage within a few degrees. This explains the transverse voltage observed in our experimental 
measurements in JJA driven near the diagonal [11] orientation. 

An analogous effect of orientational pinning has also been seen in experiments on YBCO superconductors with twin 
boundariesEj In this case, due to the correlated nature of the disorder, the direction for easy flow is the direction of the 
twins. A similar effect of horns in the parametric voltage curves are therefore observed in the direction corresponding 
to the twins. Also transport measurements when the sample is driven at an angle with respect to the twin show a 
large transverse voltage. ._. 

It is interesting to compare with the angle-dependent transverse voltage calculated for d-wave superconductors .EZI 
Also in this case, the transverse voltage vanishes only in the [10] and [11] directions. However the 9 t vs. </> curves are 
smooth in this case, since tan d t oc sin4</>.c3 This is because there is no pinning and the transverse voltage is caused 
only by the intrinsic nature of the d-wave ground state. On the other hand, the breaking of rotational symmetry 
studied here is induced by the pinning potential, and it results in non-smooth responses like "horned" parametric 
voltage curves, critical angles, transverse critical currents, etc. 

In superconductors with a square array of pinning centers, typically the pins are of circular shape and the size of 
the pins is much smaller than the distance between pinning sites.ETtj In this case, the pinning barriers that vortices 
find for motion are the same in many directions. Therefore it is possible to have orientational pinning in many of 
the square lattice symmetry directions. This explains the rich structure of a Devil's staircase observed recently in the 
simulations of Rcichhardt and Nori, where each plateau corresponds to orientational pinning in each of the several 
possible directions for orientational pinning. This interesting behavior is not possible in JJA, however, since the 
egg-carton pinning potential corresponds to the situation of square-shaped pinning centers with the pin size equal to 
the interpin distance. In this case the only possible directions for orientational pinning are the [10] and [01], as we 
have seen here. 

It is worth noting that *uany experiments in JJA in the past have been done in samples with a diagonal bias. For 
example, van Wees et al£3 have observed the existence of a transverse voltage in their measurements, which was 
unexplained. A small misalignment of their direction of the current drive could easily explain their results, since we 
have learned here that the [1 Indirection is unstable against small changes in the angle of the bias. Also Chen et a/.CJ 
have reported a transverse angle in measurements in JJA driven in the diagonal direction. In their case the effect has 
a strong component that is antisymmetric against a change in the direction of the magnetic field, which means that 
they have a Hall effect possibly due to quantum effects. However, they report that their transverse voltage had also 
a component which was even with the magnetic field (which was discounted in their computation of the Hall angle) . 
This particular spurious contribution can also be attributed to a small misalignment of the direction of the bias. From 
this we conclude that in order to study the Hall effect in JJA the most convenient choice would be a current bias in 
the [10] direction where the effect of transverse voltages at small deviations is minimum. 

When this work was upon completion, new studies of the effect of the orientation of the bias in driven square JJA 
have appeared. Fisher, Stroud and JaninEj have studied some of the effects of the direction of current in a fully 
frustrated JJA (/ = 1/2) at T = 0. In their case a transverse critical critical current and the dynamics as a function 
of I x and I y has been-described. Their results are in part complementary to our work with a single vortex (/ = 1/L ). 
Yoon, Choi and KimE3 find differences in the IV characteristics of JJA at / = when comparing the parallel current 
bias with the diagonal current bias. Their results are in agreement with our Fig. 2 results. 

In this paper we have considered the dynamics of a single vortex in a square JJA. We were able to characterize in 
detail the orientational pinning and breaking of rotational symmetry in this case. Furthermore, with the results of 
the RSJ numerical calculation we were able reproduce and interpret most of our experimental measurements for a 
quasi-diagonal bias. It remains for the future to study the behavior of a driven vortex lattice (VL) when the current is 
rotated, since the VL has also its own periodicity and symmetry directions. As we saw recently in Ref. [l^, a moving 
vortex lattice in a JJA shows different dynamical phases as a function of temperature and current. Therefore we 
expect that the characteristics of the breaking of rotational invariance, orientational pinning and transverse voltages 
will depend on the dynamical phase under consideration.!^ 
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FIG. 1. Parametric curves V y ((f>) vs. V x (</>). (a) / = 0.2, for (going outwards from the center) T = 0, T = 0.01 and T = 0.05; 
(b) I = 0.6 for T = 0.3, T = 0.5, T = 0.6, T = 1.0 and T = 1.4; (c) I = 1.2 for T = 0, T = 0.05, T = 0.1, T = 0.3 and T = 1.0. 
(d) T = 0.05 for (going outwards from the center) I = 1.0, I = 1.2, I = 1.4 and I = 1.6; (e) T = 0.2 for I = 1.0, I = 1.2, 
/ = 1.4 and I = 1.6; (f) T = 1.0 for I = 0.2, I = 0.4, I = 0.6 and I = 0.8. 

FIG. 2. Orientational pinning near <j> = 0. Voltage angle 8 V — arctan(Vg /Vx) vs. <j> for different temperatures and (a) / = 0.6 
(b) I = 1.2. 

FIG. 3. Voltage angle with respect to the [11] direction 0' v = arctan (V y — V x )/{V y + V x ) = V — n/4 vs. (f> — n/4 for different 
temperatures and (a) I = 0.6 (b) I = 1.2. 

FIG. 4. (a) Experimental results of longitudinal and transverse voltage for a current near the [11] direction, <j> = 45° ± 5 at 
T = 1.25 K. (b) Longitudinal and transverse voltage obtained numerically for 4> = 40° at T = 0.02 hIo/2ekB- (c) Idem (b) for 
an extended current range. 
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FIG. 5. (a) Experimental results of transverse angle 6 t vs. current for <fr = 45° ± 5 at different temperatures, (b) Numerical 
results of 9 t vs. / for <f> — 40°. (c) 6 t vs I at T = 0.02 for different <j> and full range of currents. Experimental temperatures are 
in Kelvin, simulation temperatures are in units of Mo/2ekB- 

FIG. 6. (a) Experimental results of transverse angle vs. T for <f> = 45° ± 5 and / = 300 fxA. (b) Numerical results of 6 t vs. 
T for = 40° and I = 0.85. 
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